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ABSTRACT 
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, The Cephei star 5 Ceti was considered one of the few monoperiodic variables 

. in the class. Despite (or perhaps because of) its apparently simple oscillation 



spectrum, it has been challenging and controversial to identify this star's pulsa- 
tion mode and constrain its physical parameters seismically. Broadband time- 
resolved photometry of 6 Ceti spanning 18.7 days with a duty cycle of about 65% 



^ ! obtained by the MOST (Microvariability & Oscillations of STars) satellite - the 

! first scientific observations ever obtained by MOST - reveals that the star is actu- 

^ ! ally multiperiodic. Besides the well-known dominant frequency of /i = 6.205886 

\0 '. d~^, we have discovered in the MOST data its first harmonic 2/i and three other 
O 
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frequencies (/s = 3.737 d-\ h = 3.673 d-^ and = 0.318 d-^), all detected 
with S/N > 4. In retrospect, /2 was also present in archival spectral line profile 
data but at lower S/N. We present seismic models whose modes match exactly 
the frequencies fi and /2. Only one model falls within the common part of the 
error boxes of the star's observed surface gravity and effective temperature from 
photometry and spectroscopy. In this model, /i is the radial {i = 0) first overtone 
and /2 is the g2 {£ = 2, m = 0) mode. This model has a mass of 10.2 ± O.2M0 
and an age of 17.9±0.3 miUion years, making S Ceti an evolved P Cephei star. If 
/2 and /a are rotationally split components of the same g2 mode, then the star's 
equatorial rotation velocity is either 27.6 km s~^ or half this value. Given its 
V sin i of about 1 km s~^, this implies we are seeing 5 Ceti nearly pole-on. 

Subject headings: stars: early- type; stars: individual (HD 16582); stars: oscilla- 
tions; stars: /3 Cephei 

1. Introduction 

Asteroseismic modeling of f3 Cephei pulsators offers an important window on the struc- 
ture and evolution of massive evolved B stars which arc precursors to corc-collapse su- 
pcrnovac. Significant progress has been made recently thanks to groundbascd single-site 
photometric campaigns lasting many years and multi-site synoptic campaigns lasting sev- 
eral months. The first method made possible the identification of 6 pulsation modes in HD 
129929 (B3V), leading to the first observational proof of non-rigid rotation inside a star other 
than the Sun (Aerts et al. 2003). The second method resulted in the identification of about 
20 frequencies in the star u Eri (B2III, Handler et al. 2004, Aerts et al. 2004, De Ridder et 
al. 2004) and a second strong case for differential internal rotation (Pamyatnykh et al. 2004) 
in a massive B star. These successes have prompted similar studies of other f] Cephei stars; 
e.g., e Oph (B2IV, Handler et al. 2005; Briquet et al. 2005) and 12 Lac (B2III, Handler et 
al. 2006). Photometry from a spacebased platform offers the advantages of both methods, 
producing long time series of much higher duty cycles than are possible from the ground. 
This paper is an indication of what is possible when such data become available for even a 
seemingly simple (3 Cephei star which has been investigated for decades with only limited 
progress. 

The bright star 5 Ceti {my — 4.07, B2IV) is one of the very few f3 Cephei stars thought 
to be a monoperiodic pulsator within that class (Stankov & Handler 2005 and references 
therein). Its variability has been investigated in a number of ground-based studies, most 
of which were based on only a few nights of data. A clear overview of these studies up to 
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1987 is not repeated here since it can be found already in Jerzykiewicz et al. (1988), who 
investigated the star's behaviour from multicolor photometry taken during seven consecutive 
nights in 1981 and one night in 1982, as well as from archival data. These authors put 
forward one oscillation frequency equal to 6.20587545 d~^ with an amplitude of ^ 12mmag 
from data assembled between 1963 and 1982. From their new data, they found night-to- 
night variations of the mean brightness and of the amplitude in the uvby filters. They 
attributed the amplitude variability to either a secondary short period with an amplitude 
below 1.6mmag or a slow drift in the data, or both. They did not find any night-to-night 
phase variations. 

Kubiak & Seggewiss (1990) collected 2 nights of simultaneous spectroscopic and photo- 
electric observations of the star and confirmed the phase lag between the radial velocity and 
light curves of 0.23 found in earher studies. 

Aerts et al. (1992) discovered large-amplitude line-profile variations in 6 Ceti, from which 
they identified the single frequency as a radial mode with a velocity amplitude of 7.4 ± 0.1 
km s~^. This identification was confirmed from available multicolor photometry by Cugier 
et al. (1994) and by Cugier & Nowak (1997). 

Finally, Daszyhska-Daszkiewicz et al. (2005) tried to put constraints on the internal 
physics of the star and to identify the radial order of the mode from combined multicolor 
photometry and radial- velocity data. They found significant differences between the data and 
their seismic models depending on whether they adopted OPAL or OP opacities, and could 
not conclude definitely if 5 Ceti was pulsating in the radial fundamental or first overtone 
mode. Neither model scenario could be made to agree with the observations. 

New space-based observations of S Ceti by the MOST satellite have now provided key 
clues to understanding this star. The first ground-based observations of 5 Ceti were obtained 
over a century ago, and over the course of the last three to four decades, there has been 
a prolonged debate about its amplitude modulation, possible multiperiodicity and period 
changes. In that time, a grand total of about 30 cycles of its dominant pulsation mode have 
been monitored. In less than three weeks, MOST was able to thoroughly sample over 70 
cycles. 



2. Observations and data reduction 

The MOST (Microvariability & Oscillations of STars) satellite (Walker, Matthews et 
al. 2003), housing a 15-cm Rumak-Maksutov optical telescope feeding a CCD photometer, 
was launched in June 2003. Its primary mission is ultraprecise rapid photometry for astero- 
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seismology. MOST was designed to monitor stars at a rate as high as 10 times per minute 
with a single-point precision of about 1-2 miUimag. For stars within the sateUite's Conti- 
nous Viewing Zone (CVZ) - in a dechnation range —18° < S < +36° - data can be collected 
nearly without interruption for up to 8 weeks. The combination of single-point precision, and 
long time coverage with high duty cycle, leads to sensitivity to rapid oscillations in Fourier 
space as low as about 1 /xmag = 1 part per million (ppm). An example of the photometric 
precision can be seen in Matthews et al. (2004). 

During the early stages of the mission when MOST was being checked out for routine 
scientific operations, the MOST Team decided to use 6 Ceti (conveniently located in the 
CVZ) as its first test target, for what was designated Commissioning Science. At the time 
of the 6 Ceti observations - the first scientific data collected by MOST - the spacecraft 
pointing had not yet been optimised, debugging of onboard software led to computer crashes 
which introduced gaps into the time series, and the downlink to Earth did not yet permit the 
nominal science data sampling rate. Hence, the photometric precision and duty cycle were 
far from what MOST was capable of achieving once commissioning was complete. However, 
at the time, this light curve represented the best combination of precision, duration and duty 
cycle ever obtained for any astronomical object other than the Sun. 

MOST monitored S Ceti for 18.68 days during 08 - 27 October 2003, with 10-s exposures 
at a sampling interval of 120 s. The data were obtained through a broadband filter designed 
expressly for the MOST instrument, with a bandpass of about 350 — 700 nm and a throughput 
approximately 3x that of a Johnson V filter. The data were reduced independently by two 
of the authors (SM and RK), yielding similar results. We present here the SM reduction, 
described below. 

The principal observing mode for MOST is known as Fabry Imaging, in which the light 
from a bright star illuminates a Fabry microlens which projects an image of the telescope 
pupil onto the CCD (see Walker, Matthews et al. 2003, and their Fig. 10). The star is 
centred in a field stop 1 arcmin in diameter, and the doughnut-shaped pupil image covers 
about 1300 pixels in a square CCD subraster of 58 pixels on a side. Seven adjacent Fabry 
microlenses sample the surrounding sky background. (One neighbouring lens is not used 
because its reading is contaminated by light from the star as pixels are transferred under the 
Target Star beam during readout.) 

Data are returned to Earth in two formats, or Science Data Streams, known as SDSl 
and SDS2. The former is processed on board so that only a small set of integrated numbers 
is sent to Earth; the latter consists of resolved Fabry images which can be fully reduced on 
the ground. Because of the limitations of the commissioning performance, the quality of the 
SDSl data (numbering about 10,000) were severely compromised and are not used for this 
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reduction. The Fabry images of the 5 Ceti SDS2 data were binned 2 x 2 to form a 29 x 29 
pixel image. For each image, all the available estimates of biases were iteratively checked 
for significant outliers and, using appropriate weights, combined to provide an average for a 
given image. 

Most of the data reduction effort is made to correct for the changing background, pri- 
marily due to orbital modulation of scattered Earthshine. Detailed aspects of this stray 
light are described by Reegen et al. (2005). Readings from the 7 adjacent Fabry microlenses 
(after 3o"-clipping of outliers) were used to derive unweighted averages of the background. 
The Primary Target Fabry field also contains about 350 pixels with only background that 
were used to calculate a third background average. This last estimate of the background 
turned out to give the best removal of stray light artifacts, and it alone was used for the 
current reduction. 

Since 5 Ceti was observed during satellite commissioning, and the pointing accuracy 
was still far from optimal, the pointing errors (with (Tx,y ~ 1 pixel) were tracked for each 
exposure, based on the ACS (Attitude Control System) telemetry. The pointing errors (and 
other possible inhomogeneities in the optics) were also estimated directly from the images 
themselves by comparing the average fluxes in the four image quadrants. Large deviations 
from quadrant to quadrant were flagged and could lead to rejection of an image. 

After bias and background subtraction, each Fabry target image was assigned a quality 
rank. This was calculated by adding 'penalties' and rejecting all images with a penalty 
count exceeding 3. Any large deviation (> 3.5(7 above average) in the X or Y pointing of the 
satellite (updated once per second) added 1 to the penalty count. A similarly large deviation 
of the (X^ + y^) pointing vector, as well as an excessive number of > 3.5(7 shifts from the 
average position in a 10-sec exposure (more than 3 such deviations in a set of 10 successive 
readings) added 1.5 each to the penalty. Any large systematic shift (> 5 pixels, either in 
X or Y) of the 10-sec average from the global average position from the entire run resulted 
in immediate rejection of the image. The rest of the quality ranking was based on dividing 
each Fabry image into four quadrants and assesing the uniformity across these quadrants. 
An iterative calculation of the percentage of outlying pixels in each quadrant (usually due 
to cosmic ray strikes or hot pixels in the CCD) sets the next component of the rank. If 
the number of rejected pixels exceeded 5% or 10%, the penalty was increased by 1 or 2, 
respectively. An image was rejected if any single quadrant contained more than 50% of the 
total number of faulty pixels (even if that number was below 10%), or any combination of 2 
quadrants was responsible for 75% of the bad pixels. 

The pre-selected pixels in each image which passed the quality control were added to 
produce an instrumental flux value. The instrument telemetry also includes the CCD focal 
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plane temperature, which was not completely regulated during the commissioning observa- 
tions. Correlations between CCD temperature and the image fluxes were corrected. The final 
stage of the reduction was high-pass filtering and Sex-clipping of the data to remove any rapid 
changes in stray light that survived the earlier quality tests. For the period range of relevance 
for 6 Ceti, this is not expected to remove any intrinsic stellar signal. As a reality check, we 
always compared the corrected fiuxes to the original raw estimates. Any large deviations 
{5m > 0.005 mag) were re-examined to avoid the danger of spurious overcorrection. 

The duty cycle of the original SDSl + SDS2 photometry was about 65%, with a stretch 
of 11 days reaching about 95%. The exclusion of the SDSl data, and the data rejection 
above, reduces the total number of data points from 3267 to 2949. The absolute duty cycle 
of this filtered light curve is only 22%, but it should be noted the reduction in duty cycle is 
primarily in each orbital cycle of 101.4 min, with no regular daily gaps as in groundbased 
photometry. The overall sampling of the S Ceti fight curve still approaches a duty cycle of 
65%. 

The total standard deviation of the 2949 data point set is 8.47 mmag, which is dominated 
by the intrinsic pulsational variability of the star. The phase diagram of the data, folded 
at the known dominant pulsation period, is presented in Fig. 1. This plot, and the light 
curve segments shown in Fig. 6, illustrate the thorough coverage and quality of the MOST 
photometry of S Ceti. 

3. Data analysis 

3.1. Discovery of the multiperiodicity of 6 Ceti 

We computed the Scargle Fourier periodogram (Scargle 1982) of the data (Fig. 2b) with 
a sampling of 10^^ d^^ and found the expected dominant peak at fi = 6.20589(8) d^^ (= 
71.827 /iHz); P = 0.16114 d = 3.867 hr. The error estimate of 8x10"^ d^^ was computed 
as (7/ = \/6 astd/7r\/iV A AT (Montgomery & O'Donohue & 1999), with N the number of 
measurements, A the amplitude, AT the total time span and tTstd the standard deviation of 
the final residuals. This value of fi is identical to within the errors, to the value reported 
by Jerzykiewicz et al. (1988). We then compared the periodogram to the spectral window 
function (Fig. 2a), which was computed from a sinusoid of the same frequency and amplitude 
as peak fi sampled at the same times as the real time series. 

A least-squares harmonic fit to the data, fixing the value of fi above, gives an amplitude 
for that peak of 11.62(3) mmag, and reveals the presence of the first harmonic 2/i with an 
amplitude of 0.72(3) mmag and a signal-to-noise of about 7.5a. While harmonics of principal 
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oscillation frequencies have been found in other (5 Cephei stars (e.g., Heynderickx et al. 1994), 
one was never before detected in 5 Ceti. This is not surprising given its amphtude is only 
6% of that of /i. 

Subtracting /i and 2/i from the data results in the periodogram shown in Fig. 2c. The 
overall variance of the time series has been reduced by 96.6%, so the residuals have a standard 
deviation of only a = 1.56 mmag. Four peaks are immediately obvious in Fig. 2c: 0.02, 2.00, 
3.73, and 14.20 d^^. The last corresponds to the known orbital frequency of the MOST 
satellite, and is an observed artifact in later MOST data, due to the modulation of scattered 
Earthshinc in the MOST focal plane with the orbital period. The 2.00 d^^ frequency is also 
almost certainly an artifact. It is due to a modulation of the stray light from the Earth 
as MOST's Sun-synchronous dusk-dawn orbit takes it over similar features of the Earth's 
albedo on a daily cycle. In this case, it may be due to the maxima in the Earth's albedo 
near both poles, not long after the Autumnal Equinox. 

The low-frequency peak near 0.02 d~^ is consistent with a long-term trend in the data, 
discussed below. The remaining peak, labeled /2 in Fig. 2c, cannot be ascribed to any known 
instrumental artifact and is almost certainly intrinsic to 5 Ceti. 

3.2. Long-term Vciriability 

In its Fabry Imaging mode, MOST is a non-differential experiment, with no comparison 
star observed in a comparable fashion in the field. However, experience with 2 years of 
MOST data has shown that the instrument is remarkably stable. Because of the previous 
observation by Jerzykiewicz et al. (1988) of a monotonic amplitude change over 7 consecutive 
nights, and the presence of a low-frequency peak in Fig. 2c, we searched for evidence of such 
amphtude changes before proceeding with frequency analysis of the MOST photometry. 

We divided the time series into subsets. Each subset was longer than the dominant 
period of 0.17 d but shorter than 0.5 d so as to mimic the previous ground-based data. Each 
subset used had at least 50 data points and no large gaps. We produced 25 such subsets and 
fitted them with /i and 2/i, allowing the amplitudes and phases to be free parameters. The 
resulting amplitudes, phases and "nightly" means of /i are plotted in Fig. 3. 

All three quantities vary beyond the error bars of the individual points, but in a com- 
plicated fashion. However, the total ranges of variability arc small: about 15% in amplitude, 
4% in phase, and about 0.4% in mean brightness. The mean brightness does show the clear- 
est evidence for a trend, which turns out to be responsible for the peak near 0.02 d~^ (P ~ 
50 d). 



-8- 



We then fitted the data residuals from Sect. 3.1 (/i and 2/i prewhitened) with a hn- 
ear trend as shown in Fig. 4. The linear fit corresponds to a brightness increase of 0.154 
mmag/day = 0.0064 mmag/hour. This is more than 2 orders of magnitude smaller than 
the trend of about 1 mmag/hour reported by Jerzykiewicz et al. (1988) in their Stromgren 
photometry of S Ceti. 

We cannot confirm whether the trend seen by MOST is stellar or instrumental (see 
discussion in Sect. 5) but we nonetheless remove it from the data for subsequent analysis. 
The detrended residuals have a standard deviation of 1.37 mmag, 0.19 mmag smaller than 
before. 



3.3. Additional oscillation frequencies in 5 Ceti 

Removing the trend from the data presented in Fig. 4 produced the periodogram shown 
in Fig. 5a. Not surprisingly, the peak at 0.02 d^^ vanished. Perhaps surprisingly, the peak 
at 14.2 d^^ also disappeared. It turned out not to be due primarily to modulated stray 
Earthshine but to the long-term trend sampled with gaps in the data at the orbital frequency 
(due to outages during the spacecraft commissioning phase). 

The peak at 2.003 d^^ does remain, consistent with a genuine stray light artifact, al- 
though its amplitude has been reduced by more than half, so the trend contributed to it as 
well. 

The frequency /2 = 3.737(2) d~^ persists as well, and its amplitude of 0.53(3) mmag - 
significant at the 5.6a level - is almost unchanged from the original data before the trend 
was removed. We note that this frequency and the dominant frequency fi have a beat period 
l/(/i — /2) = 0.405 d ~ 9.7 hr, which is well sampled by several continuous stretches in the 
MOST time series but not by the nightly stretches in the groundbased data of Jerzykiewicz 
et al. (1988). 

Prewhitening the MOST residuals by frequency /2 and the artifact at 2.003 d^^ reduces 
the standard deviation by 0.12 mmag to 1.25 mmag. From these residuals, we obtain the 
periodogram plotted in Fig. 5b. This contains two peaks with significance greater than 4(7, 
namely /g = 3.673(2) d"^ and = 0.318(2) d-^ with amplitudes of 0.39(4) (4.0(7) and 
0.43(4) mmag (4.5(7), respectively. 

Prewhitening the data by these two additional frequencies reduces the overall a by only 
0.07 mmag. A plot of the final residuals is given in Fig. 7 and their periodogram is shown in 
Fig. 5c. 
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There are peaks with amplitudes below 0.4 mmag which may help account for the 
complex amplitude and phase behaviour seen in Fig. 3. These include a peak at 3.909(5) d~^ 
at the 3.5(7 level, and one at 3.805(6) d~^ at the 2.8a level. 

3.4. Reexamination of archival data 

In hght of this frequency analysis of the MOST data, we reanalysed two high-precision 
archival data sets of S Ceti: 

1. the HIPPARCOS light curve (Ferryman et al. 1997), consisting of 72 data points 
covering 1096 days with a quasi-equidistant spacing of about 15 days; and 

2. the 60 moment variations derived from single-site high-resolution line- profile observa- 
tions by Aerts et al. (1992), spanning 7 days. 

Neither of these data sets showed any significant long-term trends. After prewhitening both 
independent data sets with frequency /i, we searched for evidence of /2 in the residuals. 
There was no sign of signal at /2 in the HIPPARCOS photometry, which is not surprising 
given its low amplitude in the MOST data and the poor samphng of the HIPPARCOS data 
for this relatively short period. However, in a periodogram of the first velocity moment 
< V > computed by Aerts et al. (1992) from their line-profile data, a peak shows up near /2, 
with an amplitude of 0.27 ±0.09 km s~^, corresponding to a significance of only 3a. Despite 
the low significance, and because of its presence in the MOST data, we conducted a second 
investigation of /2 in the archival data. 

The dominant frequency fi is present in the HIPPARCOS data with a comparable 
amplitude to that seen in the MOST data, despite the use of different custom filters for 
these observations. So we prewhitened the HIPPARCOS and MOST (detrended) lightcurves 
and the < v > moment data by /i and 2/i, and then normalised the resulting Scarglc 
periodograms to the highest peak in each. We then multiplied these normalised periodograms 
together, with the logic that if /2 (or any other common frequency) was present in more than 
one of the data sets, it would show up with improved S/N over the MOST data alone. If /2 
was absent in the other sets, it would reduce the S/N over MOST alone. In this exercise, 
all three data sets were given equal weight. 

The outcome is shown in Fig. 5c. Two peaks are immediately obvious to the eye; at 
/2 and 2.00 d~^. The latter has been reduced in signficance, consistent with it being an 
artifact in the MOST data alone, while /2 now has a significantly higher S/N of 8.6. Closer 
examination (and prewhitening) reveals the presence of a third peak at /a at the 4.6(7 level. 
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higher than in the MOST data alone. The frequency /i is not evident in the combined 
normahsed periodogram. 

These results lend circumstantial support to the presence of the frequencies /2 and /3 in 
one or both of the HIPPARCOS and line-profile data sets, with the a posteriori knowledge 
from the MOST photometry. 

3.5. Solution to the MOST light curve 

The final fit which we applied to the MOST light curve is as follows: 

6 

Ui^ a + bti + ^ Cj sm[2'K{fjti + (f)j)], (1) 

with the parameters provided in Table 1. We have retained only frequencies with amplitudes 
of significance > 4.0o", following the acceptance criterion proposed by Breger et al. (1993) 
and Kuschnig et al. (1997). 

Two segments of the fight curve comparing this solution to the data are shown in Fig. 6. 
The full line corresponds to Equation (1), while the dotted fine is a solution including only 
the linear trend and the dominant frequency /i. The residuals for the two solutions are 
shown in Fig. 7, at a magnified vertical scale. In particular, the full solution does a much 
better job near the well-populated minima and maxima of the light curve and leads to smaller 
residuals. 



4. Seismic interpretation 

4.1. Fitting fi and /2 

To compare the oscillation frequencies we have detected in 6 Ceti with those predicted 
by pulsational models of B-type stars, we explored the database described by Ausseloos et 
al. (2004). This database contains evolutionary models from the ZAMS to the turnoff with 
masses ranging from 7 to 13 M0 (in steps of O.IM©), X fixed at 0.70, Z ranging from 0.012 
to 0.030 (in steps of 0.002), and three choices for the core overshoot parameter (ctov = 0.0, 
0.1 or 0.2). For a description of the input physics of these models, and the computation of 
their oscillation frequencies, we refer to Ausseloos et al. (2004). 

Since analyses of the line-profile variations and earlier multicolor photometry indicate 
that the dominant pulsation frequency fi corresponds to a radial mode (Aerts et al. 1992; 
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Cugier et al. 1994), we have restricted our search of the model grid to those models with 
radial modes which coincide with fi . Wc have further constrained the search by forcing /2 to 
correspond to zonal modes of £2 = 0, 1 or 2, using the fitting algorithms of Ausseloos (2005). 
Modes with i > 3 are unlikely to be observed in integrated photometry due to cancellation 
effects across the stellar disk. The assumption of non-zonal modes (m 7^ 0) is considered safe 
since 5 Ceti has a low projected rotational velocity {v sin i ~ 1 km s~^; Aerts et al. 1992) so 
even a non-zero m component of mode /2 should lie close to the central component of the 
multiplet. (See our discussion in Sect. 4.2 about the inferred rotational velocity of S Ceti 
and the validity of this assumption.) 

There are no models in the database which are consistent with both fi and /2 being 
radial {i — 0) modes. We found 66 models which satisfy the criteria that /i is a radial 
mode and /2 is an £ = 1 or 2 mode; 40 models for the former, and 26 for the latter. All 
the models fitting fi and /2 simultaneously are indicated in the theoretical (log Teg, log g) 
diagram shown in Fig. 8. 

We also have independent empirical constraints on the position of 6 Ceti in Fig. 8. 
We plot the star's empirical error box derived from photometric colors averaged over the 
dominant pulsation cycle in three different systems (Walraven, Geneva and Stromgren) as 

determined by Heynderickx ct al. (1994). Wc also show the spectroscopic error box recently 
derived by Morel et al. (private communication) based on high-resolution cchelle spectra 
covering the whole pulsation cycle. Wc prefer this over previous estimates relying on a 
single spectrum during the cycle (e.g., Gics & Lambert 1994). Wc put more weight on the 
photometric error box due to its smaller extent, and the fact that it is based on three different 
photometric systems. 

It can be seen in Fig. 8 that the majority of the models matching both fi and /2 are 
too evolved (i.e., have too low gravities) for the photometric and spectroscopic error boxes. 
In fact, only five models fall within the photometric box; their physical characteristics are 
listed in Table 2. The non-radial mode corresponding to /2 is the g2 mode for all of these 
models. All these models require some amount of core overshooting. 

Which of these models is most likely to be excited? We checked the excitation rates 
with the non-adiabatic code MAD (Dupret 2001) and found that the radial first overtone of 
Model 1 in Table 2 (shown as "o" in Fig. 8) is expected to be stable due to its low metaUicity 
{Z — 0.012). All eight modes of the other four models were found to be excited. 

The high metallicity {Z = 0.028) of Model 2 (where /i is the radial fundamental mode; 
in Fig. 8) may also rule it out, since evidence points to S Ceti having a lower value of Z 
than this. Niemczura & Daszynska-Daszkiewicz (2005) recently derived [m/H] — — 0.24 ± 
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0.09 from UV iron transitions measured by lUE. Morel et al. (private communication) have 
derived the abundances of several important Z-determining elements from optical echelle 
spectra and found them to be only slightly less on average than solar values from Grevesse 
& Sauval (1998); Z = 0.02. 

We conclude therefore that the dominant mode of 6 Ceti is the radial first overtone. 

Cugicr ct al. (1994) and Cugicr & Nowak (1997) reached the same conclusion previously, 
but Daszyhska-Daszkiewicz et al. (2005) recently cast doubt on this indentification, citing 
an excitation problem. Our non-adiabatic analysis shows that this mode is unstable for 
our Models 3, 4 and 5. Of these options, wc prefer Model 3 because: (1) it is situated in 
the overlapping region of the photometric and spectroscopic error boxes, and (2) it has a 
metallicity {Z — 0.020) consistent with values derived from high-quality spectra (whereas 
the other two models appear to be too metal- rich) . The frequency /2 in this model is the g2 
i — 2 mode. 



4.2. Multiplet structure in the 5 Ceti eigenspectrum? 

What of frequency /s = 3.673 d~^? We note that it is separated from /2 by A/ = 0.064 
d~^. If /a is part of a rotationally-split multiplet, then there should be another sidelobe 
frequency at /2 + A/ = 3.801 d~^. Note that we did find a peak in the MOST data at a 
frequency of 3.805 d~^ but with a significance of only 2.8(7. 

If /2 and /a are consecutive m-values (Am = 1) of the same i — 2 quintuplet, and we 
adopt the Ledoux coefficient P-2,2 — 0.85 and radius of Model 3, we obtain a rotational 
frequency of 0.075 d~^ and an equatorial rotation velocity of 27.6 km s~^. If they are 
separated by Am = 2, then we derive half these values. 

If we combine these values with the well measured value of v sin i — 1± 1 km s~^ 
(Aerts et al. 1992), the inclination angle of 6 Ceti may be as small as i = 2° (and cannot be 
larger than i — 8°). Hence we must be observing the star nearly pole-on. This conclusion is 
independent of the choice of models in Table 2. 

The pole of a star corresponds to an angle of complete cancellation for an £ = 2 sectoral 
and tesseral mode (Chadid et al. 1999, their Table A. 1). It is therefore most hkely that /2 
corresponds to the central peak of the quintuplet, as assumed in the modeling. Even if it 
were the outermost component of the quintuplet, the central peak would differ only by ~ 0.13 
d~^ from /2. Such a frequency shift is not large enough to affect our model identifications 
which assumed a low rotation rate and hence, closely spaced multiplet structure. 
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We also called attention in Sect. 3.3 to another frequency, at 3.909 d"^, with a signifi- 
cance of 3.5(7. This frequency would not fit into an equidistant quintuplet structure which 
includes /2 and f^. We do note that Model 3 has a. £ — 3 mode predicted to be excited 
whose frequency is close to 3.909 d~^. 

Finally, we have no obvious explanation for frequency /i = 0.318 d"^, but we do point 
out that /4 ~ 5(/2 — /a) ~ 3(3.909 — 3.805) d~^, within the errors. Perhaps it is related to 
beating between modes. 

5. Summary 

MOST photometry of the P Cephei star S Ceti - until now regarded as a prototypical 
example of a monoperiodic radial oscillator within the class - reveals the presence of at least 
two additional oscillations consistent with non-radial modes, as well as the first harmonic of 
the dominant radial mode. We have compared the new frequency spectrum with pulsation 
models constrained by stellar parameters based on photometric colors and spectroscopic 
analysis. Wc conclude that the dominant mode in 6 Ceti is due to the radial first overtone, 
and that the next strongest mode is the g2 i = 2 mode. We investigated multiplet structure 
associated with the latter mode to constrain the rotational velocity of the star, and showed 
that 5 Ceti is seen nearly pole-on. 

We also find a shallow linear brightness increase in the star, at a rate of about 0.0064 
mmag/hour. Jerzykiewicz et al. (1988) also found a drift in their ubvy photometry of 5 Ceti, 
but at a rate of about 1 mmag/hour. If both trends were truly intrinsic to the star, then the 
rate of brightness change would be strongly variable as a function of epoch. Such variations 
would not be surprising for an evolved star like 5 Ceti undergoing small instabilities on its 
way to the end of the central hydrogen burning phase. However, it is still possible that 
the trend seen by Jerzykiewicz et al. (1988) was an uncorrected effect of extinction in the 
Earth's atmosphere. MOST does not suffer such extinction effects in orbit, but as a non- 
differential experiment, it is impossible to exclude an instrumental origin for the gradual 
trend it measured. 

We note that Jerzykiewicz et al. (1988) also observed amplitude modulation in their 
night-to-night photometry which they suggested might be due to a second short-period 
variation with an amplitude below 1.6 mmag. The detection by MOST of frequency /2, with 
an amplitude of about 0.5 mmag is consistent with that explanation. 

The best model fit indicates that 6 Ceti has a mass of M = 10.2 ± 0.2Mq, an age 
of 17.9±0.3 million years and core overshooting (ojov = 0.20 ± 0.05). It is only the third 
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(3 Cephei star for which the core overshooting parameter could be determined, along with 
HD 129929 (ttov = 0.10 ± 0.05, Aerts et al. 2003) and p Eridani (ctov < 0.13, Pamyatnykh 
et al. 2004). The determination of the core overshooting parameter results in an accurate 
seismic mass estimate, which is for all three stars an improvement by an order of magnitude 
over previous photometric or spectroscopic mass estimates. 

These results illustrate the power of nearly continuous ultra-precise photometry in un- 
derstanding the structure and evolution of f] Cephei stars. We can anticipate additional 
results from the space-based observatories WIRE (Briintt et al., in preparation), MOST, 
COROT and Kepler in the coming years. 
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from the Natural Sciences & Engineering Research Council (NSERC) Canada. RK's work 
is supported by the Canadian Space Agency. 
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Phase 



Fig. 1. — Phase diagram of the MOST data of 5 Ceti folded at the known dominant pulsation 
period of 0.161138 days 
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Fig. 2. — Periodograms of the MOST photometry of 5 Ceti. Panel a: window function shifted 
and scaled to the main peak; panel b: periodogram of the data; panel c: periodogram after 
prewhitening with /i and 2/i. 



18 




HJD - 2451300. 




HJD - 2451300. 



Fig. 3. — Amplitude of /i, phase of /i and mean with respect to the overall mean, for 25 
datastrings observed by MOST. 
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Fig. 4. — Brightness increase of 6 Ceti in the residuals after prewhitening the dominant 
frequency /i and its harmonic from the MOST hghtcurve. 



-20- 



□ 
E 

T3 

"q. 
E 
< 



cri 
o 
E 

"O 

"5. 
E 
< 




in 
d 




10 15 
Frequency (c d~^) 



20 



Fig. 5. — Periodograms of S Ceti. Panel a: after prewhitening with fi and 2/i and detrending 
according to Fig. 4. Panel b: after subsequent prewhitening with /2 and 2.003 d~^. Panel c: 
product of normalised amplitude spectra of the MOST and HIPPARCOS photometry and 
the first velocity moment < v > taken from Aerts et al. (1992). 
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Fig. 6. — Comparison between the MOST data and (1) the final fit given in Eq. (1) (full 
line), (2) a fit including only the dominant frequency and the trend (dotted line) for a few 
selected segments. 
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Fig. 7. — Residuals of the entire MOST data set after prewhitening by the final fit given 
in Eq. (1) (upper panel) and by only the dominant frequency /i and the linear trend in the 
data (lower panel). The standard deviations of these residuals are 1.18 and 1.47 mmag, 
respectively. 
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Fig. 8. — Position of stellar models in the (log Teg, log 51) diagram whose zonal modes fit /i 
and /2 simultaneously. The symbol convention is as follows: fi radial first overtone and /2 
an i — l,m. — g mode (o), fi radial second overtone and f2 an i — 1,111 — g mode (□), 
/i radial fundamental and f2 an i — 2,171 — g mode (+), /i radial first overtone and /2 
an i = 2,m = g mode (x). The photometric (full lines) and spectroscopic (dotted lines) 
observational error boxes of S Ceti are also indicated. 
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Table 1. Final lightcurve solution of S Ceti, according to the eight terms given in Eq. (1). 
The reference epoch for the phases is the time of the first measurement 

(HJD 2451377.140428). 



fj (cd-^) 


(/.Hz) 


Cj (mmag) 




S/N 


/i = 6.20589(8) 


71.8274(9) 


11.62(3) 


0.4958(4) 


121.4 


2/i 


2/i 


0.72(3) 


0.619(7) 


7.5 


2.003(1) 


23.18(1) 


0.65(3) 


0.292(7) 


6.8 


/2 = 3.737(2) 


43.25(2) 


0.53(3) 


0.454(9) 


5.5 


/3 = 3.673(2) 


42.51(2) 


0.39(4) 


0.82(1) 


4.0 


U = 0.318(2) 


3.68(2) 


0.43(4) 


0.55(2) 


4.5 


a = 922.89(5) mmag 



b = —0.154(5) mmag/day 
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Table 2. Physical parameters of the stellar models within the photometric error box of 
5 Ceti shown in Fig. 8. The core overshoot parameter is expressed in units of the local 

pressure scale height. 





log L/Lq 


log T^ff 




M(M0) 


R{Re) 




X 


Z 


Q^ov 


age(yr) 


l:o 


3.977 


4.332 


3.714 


9.43 


7.06 


0.088 


0.70 


0.012 


0.0 


19.7x10^ 


2:+ 


4.014 


4.334 


3.727 


10.36 


7.29 


0.252 


0.70 


0.028 


0.1 


16.9x10^ 


3:x 


4.064 


4.347 


3.725 


10.23 


7.27 


0.249 


0.70 


0.020 


0.2 


17.9x10^ 


4:x 


4.026 


4.339 


3.721 


10.04 


7.23 


0.255 


0.70 


0.022 


0.2 


18.6x10^ 


5:x 


3.992 


4.331 


3.718 


9.88 


7.20 


0.260 


0.70 


0.024 


0.2 


19.3x10^ 



